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The  efficacy  of  dielectric  barrier  discharge  (DBD)  plasmas  driven  by  high  voltage  (~15 
kV)  repetitive  nanosecond  pulses  (~100  ns  FWHM)  for  flow  separation  control  is 
investigated  experimentally  on  an  airfoil  leading  edge  up  to  Re=lxl06  (62  m/s).  Unlike  AC- 
DBDs,  the  nanosecond  pulse  driven  DBD  plasma  actuator  transfers  very  little  momentum  to 
the  neutral  air,  but  generates  compression  waves  similar  to  localized  arc  filament  plasma 
actuators.  A  complex  pattern  of  quasi-planar  and  spherical  compression  waves  is  observed 
in  still  air.  Measurements  suggest  that  some  of  these  compression  waves  are  generated  by 
discharge  filaments  that  remain  fairly  reproducible  pulse-to-pulse.  The  device  performs  as 
an  active  trip  at  high  Re  pre-stall  angles  of  attack  and  provides  perturbations  that  generate 
coherent  spanwise  vortices  at  post-stall.  These  coherent  structures  entrain  freestream 
momentum  thereby  reattaching  the  normally  separated  flow  to  the  suction  surface  of  the 
airfoil.  Coherent  structures  are  identified  at  all  tested  frequencies,  but  values  of  Fc+= 4-6  are 
most  effective  for  control.  Such  devices  which  are  believed  to  function  through  thermal 
effects  could  be  an  alternative  to  AC-DBD  plasmas  that  rely  on  momentum  addition. 
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Nomenclature 


a 

- 

local  speed  of  sound 

c 

= 

model  chord,  25.4  cm 

d 

= 

plasma  filament  diameter 

f 

= 

frequency 

%sep 

= 

separation  length  scale 

u 

= 

instantaneous  streamwise  velocity 

u 

= 

ensemble-averaged  streamwise  velocity 

u 

= 

phase-averaged  streamwise  velocity 

f/oo 

= 

freestream  velocity 

V 

= 

instantaneous  streamwise  velocity 

V 

= 

ensemble-averaged  normal  velocity 

V 

= 

phase-averaged  normal  velocity 

X 

= 

streamwise  coordinate 

y 

= 

vertical  coordinate 

a 

= 

angle  of  attack  in  degrees 

P 

= 

air  density 

^ rise 

= 

pulse  rise  time 

^ acoustic 

= 

acoustic  time  scale 

X 

= 

plasma  wavelength 

^ci 

= 

swirling  strength 

4,* 

= 

normalized  swirling  strength,  Xcic/U 

V 

= 

kinematic  viscosity 

cD 

= 

drag  coefficient 

CL 

= 

sectional  lift  coefficient 

cP 

= 

fluctuating  pressure  coefficient 

cP 

— 

pressure  coefficient 

2 

American  Institute  of  Aeronautics  and  Astronautics 

F  + 

1  c 

Re 
0 
Q 
Q * 


=  reduced  frequency  based  on  model  chord,  fc/U^ 
=  chord  based  Reynolds  number,  U^c/v 
=  phase 

=  spanwise  component  of  vorticity 
=  normalized  spanwise  vorticity,  Qc/U ^ 


I.  Introduction 

Flow  over  solid  bodies  is  retarded  near  the  surface  due  to  viscous  effects  which  give  rise  to  a  boundary  layer. 
Often  the  near  wall  fluid  must  also  negotiate  an  adverse  pressure  gradient  and/or  geometric  aberration.  Flow 
separation  occurs  when  the  momentum  in  the  boundary  layer  is  insufficient  to  overcome  these  impediments.  This  is 
marked  by  a  thickening  of  the  rotational  region  of  fluid  near  the  surface  and  an  increase  in  the  wall  normal  velocity 
component.  Boundary  layer  separation  is  synonymous  with  substantial  energy  losses  and  minimizing  its  occurrence 
is  one  of  the  primary  design  objectives  for  aerodynamic  and  hydrodynamic  systems. 

A  common  platform  for  investigating  separation  control  is  the  airfoil  at  high  angle  of  attack.  This  system 
produces  maximum  lift  near  the  angle  at  which  separation  occurs  from  the  leading  edge.  Slight  excursions  to  higher 
angles  can  result  in  unrecoverable  stall.  Transport  aircraft  often  combat  this  problem  using  passive  control  devices  in 
the  form  of  segmented  leading  edge  slats  that  deploy  during  takeoff  and  landing  to  promote  mixing  between  high 
momentum  fluid  on  the  pressure  surface  with  low  momentum  fluid  on  the  suction  surface.  These  passive  geometric 
modifications  are  effective  for  augmenting  lift,  but  create  significant  increases  in  mechanical  complexity, 
manufacturing  cost,  weight  and  parasitic  drag  even  when  stowed  during  cruise. 

Active  flow  control  in  the  form  of  localized  periodic  excitation  has  been  used  to  impart  global  changes  to 
separated  flow  fields.1  This  is  often  accomplished  by  adding  spanwise  slots  to  the  aerodynamic  surface  near 
separation  locations.  Periodic  excitation,  usually  in  the  form  of  momentum,  is  delivered  through  the  slots  to 
manipulate  natural  flow  instabilities  in  an  effort  to  prevent  or  reattach  separated  flows.  The  optimal  dimensionless 
frequency,  F^=fxsep/U0 0,  for  controlling  separation  using  periodic  excitation  is  on  the  order  of  unity  where/  xsep  and 
Uoo  are  the  forcing  frequency,  separation  length  scale  and  freestream  velocity  respectively.2  Periodic  excitation  can 
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be  accomplished  through  a  variety  of  mechanisms,  but  piezoelectric,  electromagnetic  or  electrostatic  devices  are 
often  employed.  In  ah  of  these  cases,  an  electromechanical  driver  creates  the  oscillatory  flow  used  for  excitation. 
These  devices  are  controlled  through  electrical  signals  and,  compared  to  passive  control,  offer  a  significant  reduction 
in  weight,  mechanical  complexity  and  parasitic  drag.  Unfortunately,  they  possess  limited  bandwidth  and  are  subject 
to  mechanical  failure  because  the  electromechanical  driver  is  usually  operated  at  resonance  to  produce  the  large 
amplitude  perturbations  necessary  at  practical  flight  speeds.  Even  when  operated  in  this  fashion,  it  is  difficult  to 
meet  amplitude  requirements  especially  for  cruise  conditions  since  the  required  actuator  momentum  scales  roughly 
with  dynamic  pressure.1 

Flow  control  with  plasma  actuation  is  appealing  because  these  devices  are  surface  mounted,  lack  mechanical 
parts  and  require  relatively  low  power.  Dielectric  Barrier  Discharge  (DBD)  plasma  actuators  driven  by  AC 
waveforms  (AC-DBD)  are  the  most  popular.3  They  have  been  demonstrated  for  airfoil  leading  edge  separation 
control  at  freestream  velocities  up  to  60  m/s  (Re~  106).4  The  control  mechanism  for  AC-DBD  plasma  actuators  arises 
from  an  electrohydrodynamic  (EHD)  effect  (i.e.  Coulomb  force  interaction).  The  signature  of  this  mechanism  is 
generation  of  a  low  speed  (usually  <10  m/s)  near  wall  jet  in  quiescent  air  due  to  interaction  between  the  charged 
species  in  the  plasma  and  neutral  particles  near  the  surface.  The  technology  continues  to  improve  through 
optimization  of  a  single  actuator,  the  use  of  multiple  actuators/electrodes  and  sliding  discharges  that  rely  on 
additional  DC  potential.5"8  The  common  goal  in  these  cases  is  to  increase  the  momentum  production  of  the  device. 

Early  reports  suggest  DBD  plasma  actuators  driven  by  a  different  type  of  waveform  could  be  an  alternative.  The 
construction  of  the  device  is  analogous  to  the  AC-DBD,  but  the  discharge  is  driven  by  repetitive  nanosecond 
duration  pulses  (FWHM-5-100  ns)).  We  will  refer  to  this  type  of  actuation  as  a  repetitive  nanosecond  pulse 
discharge  or  “NS-DBD”  throughout  the  paper.  DBD  plasma  created  using  these  waveforms  has  shown  control 
authority  for  leading  edge  airfoil  separation  control  up  to  Mach  0.8 5, 9  but  many  questions  about  the  underlying  flow 
physics  remain  unanswered.  It  should  be  noted  that  such  waveforms  have  also  been  used  in  combination  with  AC 
and  DC  bias  voltage,  but  to  the  authors  knowledge  these  have  not  been  tested  in  flow  control  studies.10,11 

The  mechanism  by  which  NS-DBD  plasmas  produce  control  is  believed  to  be  substantially  different  from  the 
AC  discharge.  The  NS-DBD  produces  very  low  velocity  in  the  neutral  species.  Instead,  a  thermal  effect  (i.e.  Joule 
heating)  that  generates  a  local  compression  wave  is  believed  to  provide  the  control  mechanism.  Indeed,  compression 
waves  generated  in  quiescent  air  by  localized  heating  of  the  near  surface  gas  layer  have  been  detected  in  both 
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modeling  and  experiment.9,12  This  scenario  is  well-established  for  localized  arc  filament  plasma  actuators  (LAFPAs) 
that  are  effective  for  controlling  high  Reynolds  number  and  high-speed  (subsonic/supersonic,  cold/hot)  jets  in  both 
experiments13"15  and  computations.16  Like  many  flow  control  studies,  the  main  idea  of  the  LAFPA  approach  is  to 
force  the  flow  at  a  frequency  approaching  one  of  the  flow  instabilities,  thereby  triggering  subsequent  growth  of 
seeded  perturbation  and  their  rollup  into  flow  structures.  In  high  speed  flow  this  requires  both  high  amplitude  and 
high  bandwidth  to  realize  performance  goals  across  various  flight  regimes.  Studies  of  flow  control  using  LAFPAs  in 
atmospheric  pressure  cold  and  hot  jets  (M=0. 9-2.0)  demonstrate  significant  localized  heating  and  repetitive 
compression  wave  formation  which  excite  flow  instabilities  that  alter  dynamics  of  large-scale  coherent  structures 
thereby  achieving  noise  reduction  and  mixing  enhancement.  In  fact,  these  are  the  only  devices  that  have  shown  this 
capability  at  high  speeds.  This  has  been  accomplished  with  relatively  low  actuator  power,  (-10  W  per  actuator),  at 
forcing  frequencies  near  the  jet  column  (preferred  mode)  instability  frequency.  Thus,  the  main  benefit  for  flow 
control  with  thermal  perturbations  seems  to  be  their  use  at  high  speeds  where  it  is  difficult  for  momentum  based 
devices  to  reliably  or  efficiently  produce  sufficient  momentum  for  control  authority.  It  should  be  noted  that 
compared  to  the  EHD  effect,  the  effect  of  Joule  heating  in  surface  DBD  actuators  driven  by  AC  voltage  waveforms 
appears  insignificant.17'19 

This  paper  describes  the  development  and  application  of  nanosecond  pulse  driven  DBD  plasma  actuators  for 
controlling  flow  separation  from  an  airfoil  leading  edge.  The  NS-DBD  plasma  actuator  is  characterized  using 
electrical  measurements  and  optical  techniques  in  quiescent  air.  Control  authority  is  demonstrated  on  an  airfoil 
leading  edge  up  to  Re=  lxl 06  (62  m/s)  using  static  pressure  and  PIV  measurements.  Results  show  the  NS-DBD  can 
perform  as  an  active  boundary  layer  trip  or  generate  coherent  spanwise  vortices  that  are  widely  associated  with 
controlling  separation  in  a  variety  of  flow  systems.1  The  employed  control  mechanism  (active  trip  or  coherent 
vortex)  is  dependent  on  the  Reynolds  number  and  angle  of  attack.  The  generation  of  large-scale  structures  via 
thermal  perturbations  is  in  agreement  with  experiments  using  LAFPAs  on  hot  and  cold  high  speed  high  Re  jets.  This 
suggests  that  the  mechanisms  of  flow  control  for  both  types  of  actuators  may  be  similar. 

II.  Experimental  Facilities  and  Techniques 

A.  NS-DBD  Plasma  Actuator 
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The  DBD  plasma  actuators  in  this  work  are  composed  of  two  copper  tape  electrodes  separated  by  Kapton  tape 
dielectric  arranged  in  an  asymmetric  fashion  (Figure  1).  There  is  a  slight  overlap  between  the  electrodes  to 
encourage  uniform  plasma  generation.  The  covered  ground  electrode  is  12.7  mm  (lA  inch)  wide  and  the  exposed 
high  voltage  electrode  is  6.35  mm  (%  in)  wide.  Both  electrodes  have  thickness  of  0.09  mm  (3.5  mil).  The  dielectric 
barrier  is  composed  of  3  layers  of  Kapton  tape.  Each  layer  has  thickness  of  0.09  mm  (3.5  mil)  and  dielectric  strength 
of  10  kV.  Each  layer  of  Kapton  tape  has  a  0.04  mm  (1.5  mil)  layer  of  silicone  adhesive  such  that  the  actual  Kapton 
thickness  for  each  tape  layer  is  only  0.05  mm  (2  mil).  The  total  thicknesses  of  the  dielectric  and  the  device  as  a 
whole  are  0.27  mm  (10.5  mil)  and  0.44  mm  (17.5  mil)  respectively.  The  actuator  construction  has  been  determined 
based  on  previous  studies  of  AC-DBD  plasma  and  has  not  been  optimized  for  NS-DBD  purposes.20,21  The  actuator  is 
directly  applied  to  the  substrate  in  all  cases. 


High  voltage  electrode 

6.35mm  wide  0.1mm  thick\ 

Three  layers  of  0.1  mm  Kapton®  tape 

i _ !_ 

^ _ 

(TO  L  J 

- 3 — u  Ground  electrode 

12.7mm  wide 

s.  0.1  mm  th  ick 

r< - T\ 

'Poise generator  V  ' 

P6015A  ^ 

High  voltage  probe 

< 

< 

^  0.2Q  Current  shunt  resistor 
>  with  49. 8Q  output  resistance 

— 

Figure  1.  DBD  Plasma  actuator  schematic  with  connections  for  electrical  measurements 


Two  surface  DBD  plasma  setups  have  been  used  for  actuator  characterization.  The  first  setup,  shown  in  Figure 
1,  is  mounted  on  a  flat  dielectric  surface  and  used  for  nanosecond  discharge  pulse  energy  and  induced  flow 
measurements.  The  second  setup,  used  for  schlieren  imaging,  employed  an  actuator  with  the  same  electrode 
geometry  mounted  on  a  5  cm  diameter  cylinder.  In  these  measurements,  the  actuator  length  (into  the  page)  is  fixed  at 
30  cm.  The  cylinder  is  chosen  to  simplify  alignment  of  the  schlieren  diagnostic.  Figure  2  shows  a  photograph  of  the 
second  actuator,  with  the  repetitive  nanosecond  pulse  discharge  in  operation.  In  all  cases,  the  repetitive  nanosecond 
pulse  voltage  waveform  for  the  surface  DBD  discharge  is  generated  by  high-voltage  magnetic  compression  type 
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pulsed  power  supplies,  custom  designed  and  built  at  The  Ohio  State  University.  Briefly,  these  power  supplies 
generate  high  voltage  pulses  (peak  voltage  up  to  20  kV,  pulse  duration  50-100  nsec  FWHM),  at  a  continuous  pulse 
repetition  rate  of  up  to  3  kHz,  and  pulse  energy  of  up  to  100  mJ/pulse.  Note  that  there  is  no  lower  bound  on  the 
repetition  frequency  and  no  low  frequency  burst  or  amplitude  modulation  is  performed  in  this  work.  Both  the  peak 
voltage  and  the  pulse  energy  coupled  to  the  load  are  strongly  load-dependent.  The  pulse  waveform  is  transmitted  to 
the  DBD  actuator  load  via  a  collinear  pair  transmission  line  up  to  1.5  m  long. 


Figure  2.  Photograph  of  a  NS-DBD  actuator  mounted  on  a  cylinder  in  operation 
B.  Plasma  Diagnostics 

Pulse  voltage  and  current  have  been  measured  using  a  Tektronix  P6015  high  voltage  probe,  a  custom  made 
current  shunt  probe,  and  a  Tektronix  TDS2024  oscilloscope.  Both  are  measured  at  the  end  of  the  transmission  line, 
where  the  transmission  line  cables  are  connected  to  the  load.  The  current  probe  uses  five  resistors  of  1  Q  each,  with 
the  total  resistance  of  0.2  f},  arranged  radially  to  minimize  residual  inductance.  The  impedance  of  the  current  probe 
is  set  at  50  ^  by  adding  a  49.8  ^  resistor.  A  schematic  of  the  setup  is  shown  in  Figure  1.  For  termination  of  the 
current  probe  signal,  a  50  Q  feed-through  terminator  resistor  is  used  at  the  oscilloscope.  This  reduced  multiple 
reflections  of  the  current  probe  signal  from  the  oscilloscope  and  from  the  current  probe.  The  stray  phase  shift 
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between  the  high  voltage  probe  signal  and  the  current  probe  signal  is  determined  by  reducing  to  zero  the  integral  of 
the  product  of  voltage  and  current  at  the  conditions  when  the  pulse  voltage  is  low  (a  few  hundred  Volts)  (i.e.  when 
breakdown  between  the  actuator  electrodes  is  not  produced).22 

Compression  waves  generated  by  the  nanosecond  pulse  discharge  in  the  DBD  actuator  are  detected  using  a 
custom-build  phase-locked  schlieren  system.  The  light  source  is  a  5  W  green  light  emitting  diode  (LED)  controlled 
by  a  custom  made  drive  circuit,  which  can  operate  both  in  continuous  emission  and  pulsed  mode.  To  enhance  the 
sensitivity,  a  short  focal  length  plano-convex  lens  is  used  to  focus  the  LED  light  at  the  knife  edge.  The  schlieren 
images  are  acquired  by  a  Thorlab  DCC1545  CMOS  camera.  To  discriminate  stray  UV/visible  emission  from  the 
plasma,  a  green  dichroic  filter  (90%  band  pass  transmission  at  510-570  nm  for  0°  angle  of  incidence)  is  placed  in 
front  of  the  camera.  The  maximum  transmission  wavelength  is  varied  by  rotating  the  filter.  The  LED  pulse  duration 
is  set  at  300  ns  which  gives  the  required  MHz  time  resolution  and  high  signal  to  noise  ratio  necessary  for  resolving 
the  generated  compression  waves.  The  LED  is  operated  in  a  repetitively  pulsed  mode  at  the  same  pulse  repetition 
rate  as  the  nanosecond  pulse  discharge  in  the  DBD  actuator  (1  kHz).  The  LED  drive  circuit  is  synchronized  with  the 
pulsed  plasma  generator  with  a  certain  delay  time  and  phase-locked  schlieren  images  are  averaged  over  150  to  220 
discharge  pulses.  A  reference  time  is  defined  as  the  moment  when  the  discharge  current  reaches  50%  of  its  peak 
value.  The  delay  time  between  this  reference  and  the  LED  emission  pulse  is  approximately  3  psec.  The  pulsed 
plasma  generator  jitter  is  approximately  50-100  nsec  and  with  the  overall  timing  uncertainty  of  ±200  nsec  the 
uncertainty  of  the  wave  front  position  becomes  0.07  mm. 

C.  Wind  Tunnel  Facility  and  Airfoil  Model 

A  recirculating  wind  tunnel  with  an  optically  accessible  61x61x122  cm3  (2x2x4  ft3)  closed  test  section 
serves  as  the  test  bed.  Test  sections  walls  are  constructed  of  25.4  mm  (1  in)  thick  super  abrasion  resistant  acrylic. 
Each  side  wall  is  fitted  with  a  30.5  cm  (12  in)  diameter  port  that  is  located  30.5  cm  (12  in)  from  the  test  section  floor 
and  61  cm  (24  in)  downstream  of  the  test  section  entrance.  Air  flow  in  the  tunnel  is  continuously  variable  from  3-90 
m/s  (10-300  ft/s).  Flow  conditioning  upstream  of  the  test  section  includes  a  hexagonal  cell  aluminum  honeycomb 
while  high  porosity  stainless  steel  screens  are  mounted  downstream  of  the  test  section  as  a  safety  catch.  Four  high 
efficiency  turning  cascades  fabricated  of  galvanized  steel  are  installed  in  each  of  the  four  tunnel  elbows.  This 
assembly  results  in  freestream  turbulence  levels  of  -0.25%  near  the  test  section  inlet.  The  tunnel  is  also  equipped 
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with  a  commercial  aluminum  fm/copper  tube,  double  row  heat  exchanger  with  set  point  controller  and  electronic 
modulating  valve.  This  arrangement  allows  the  tunnel  freestream  operating  temperature  to  be  maintained  within  1°C 
from  the  ambient  when  supplied  a  sufficient  source  of  cooling  water  (max  189  1pm  (50  gpm)). 

The  airfoil  is  a  simplified  high-lift  version  of  the  NASA  Energy  Efficient  Transport  (EET)  with  chord  length  of 
25.4  cm  (10  in)  and  span  of  61  cm  (2  ft).  It  was  originally  designed  for  studies  of  separation  control  over  a  deflected 
flap  at  low  angle  of  attack.20,21  The  current  work  examines  separation  from  the  leading  edge  at  high  angle  of  attack 
with  no  flap  deflection.  The  model  is  constructed  of  a  nylon  compound  (Duraform  GF)  and  has  been  fabricated 
using  selective  laser  sintering  (SLS)  technology.  Angle  of  attack  settings  are  done  manually  using  a  protractor  on  the 
side  wall.  Forty  five  static  pressure  taps  are  staggered  around  the  test  section  centerline.  The  airfoil  is  mounted 
across  the  span  of  the  tunnel  and  supported  by  12  inch  diameter  wall  plugs  on  either  side.  Maximum  blockage  in  the 
current  study  is  -12%  and  occurs  at  16°  angle  of  attack.  No  blockage  or  wall  corrections  have  been  employed.  While 
it  is  quite  challenging,  future  studies  should  investigate  the  effect  of  blockage  on  active  flow  control  in  general.  To 
the  author’s  knowledge  this  has  not  been  addressed  in  the  literature.  Significant  studies  on  active  flow  control  with 
zero  net  mass  flux  actuators  have  been  completed  for  a  more  complex  version  of  this  airfoil.23"28 

D.  Flow  Diagnostics 

Measurements  of  static  pressure  from  taps  on  the  model  surface  are  acquired  using  Scanivalve  digital  pressure 
sensor  arrays  (DSA-3217).  Subsequent  values  of  sectional  CP  and  CL  are  averaged  over  50  samples  acquired  at  1  Hz. 
A  flush  mounted  Kulite  pressure  transducer  (XCQ-080-25A)  in  the  model  at  x/c=0.40  is  powered  using  an  in-house 
constructed  signal  conditioner  that  amplifies  output  by  1000.  A  separate  Kemo  unit  is  employed  for  low-pass 
filtering  at  25  kHz.  The  resulting  pressure  traces  are  sampled  simultaneously  at  50  kHz  using  a  National  Instruments 
PCI-6143  data  acquisition  board.  Average  spectra  are  calculated  from  50  blocks  of  8192  samples  which  results  in  a 
frequency  resolution  of  approximately  6  Hz.  A  Hanning  window  function  is  applied  to  each  block. 

Two-component  particle  image  velocimetry  (PIV)  is  used  to  obtain  quantitative  measurements  of  the  velocity 
fields  over  the  airfoil  and  the  DBD  actuator  on  a  flat  plate  in  still  air.  Images  are  acquired  and  processed  using  a 
LaVision  PIV  system.  Nominally  submicron  olive  oil  seed  particles  are  introduced  upstream  of  the  test  section 
contraction  using  a  6-jet  atomizer.  A  dual-head  Spectra  Physics  PIV-400  Nd:YAG  laser  is  used  in  conjunction  with 
spherical  and  cylindrical  lenses  to  form  a  thin  light  sheet  that  allows  PIV  measurements.  The  time  separation 
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between  laser  pulses  used  for  particle  scattering  is  set  according  to  the  flow  velocity,  camera  magnification  and 
correlation  window  size.  Two  images  corresponding  to  the  pulses  from  each  laser  head  are  acquired  by  a  LaVision 
14  bit  2048  by  2048  pixel  Imager  Pro-X  CCD  camera  equipped  with  a  Nikon  Nikkor  50  mm  f/1.2  lens.  Flat  plate 
PIV  data  are  acquired  using  an  additional  2x  teleconvertor  to  increase  spatial  resolution.  For  each  image  pair, 
subregions  are  cross-correlated  using  decreasing  window  size  (642-322  pixel2)  multi-pass  processing  with  50% 
overlap.  The  resulting  velocity  fields  are  post-processed  to  remove  spurious  vectors  using  an  allowable  vector  range 
and  median  filter.  Removed  vectors  are  replaced  using  an  interpolation  scheme  based  on  the  average  of  neighboring 
vectors.  A  3x3  Gaussian  smoothing  filter  is  also  applied  to  the  calculated  velocity  fields.  The  PIV  data  are  sampled 
at  10  Hz.  Ensemble-averaged  statistics  for  the  airfoil  and  actuator  characterization  are  calculated  from  750  and  500 
instantaneous  velocity  fields,  respectively.  Convergence  calculations  show  these  values  are  more  than  sufficient  for 
the  data  presented  here. 

Phase-locked  PIV  data  are  acquired  using  the  programmable  timing  unit  of  the  LaVision  system.  The  acquisition 
is  synced  with  the  frequency  of  the  actuation  signal.  Velocity  fields  at  four  phases  of  the  actuator  frequency  are 
investigated  by  stepping  through  the  forcing  period  using  time  delays.  The  resulting  phase-locked  data  sets  are 
averaged  over  100  images  at  each  phase  which  convergence  calculations  show  are  sufficient  for  resolving  the 
primary  features  of  the  average  flow  fields  (i.e.  coherent  structures).  Phase-locked  PIV  data  is  acquired  at  5  Hz.  The 
spatial  resolutions  of  PIV  data  for  the  airfoil  and  flat  plate  data  sets  are  approximately  2.4  and  0.2  mm,  respectively. 
Near  surface  measurements  for  the  flat  plate  data  sets  are  obtained  within  0.3  mm  of  the  substrate. 

Full-scale  accuracy  for  measurements  of  instantaneous  velocity  is  1.7%  for  actuator  characterization  and  0.9%  in 
the  case  of  the  airfoil.  This  is  calculated  by  assuming  negligible  laser  timing  errors  and  a  correlation  peak  estimation 
error  of  0.1  pixels.  The  error  is  smaller  in  the  airfoil  case  due  to  better  choices  of  time  separation  between  laser 
pulses  due  to  the  a  priori  knowledge  of  maximum  velocities  in  the  freestream  flow.  The  uncertainties  in 
measurements  of  mean  velocity  are  calculated  using  95%  confidence  intervals.29  This  results  in  relative  error  of  4- 
6%  for  actuator  characterization  and  2-3%  in  the  airfoil  case.  The  value  for  the  actuator  in  still  air  case  is  estimated 
using  the  average  standard  deviation  and  average  velocity  in  the  near  wall  profiles  which  gives  a  more  representative 
estimate  than  the  maximum  velocity.  The  uncertainty  in  the  airfoil  case  is  estimated  using  the  maximum  standard 
deviation  in  the  turbulent  wake  and  the  freestream  velocity.  A  range  of  values  that  encompasses  all  the  flow  fields 
surveyed  is  provided  in  the  interest  of  brevity.  The  95%  confidence  interval  for  CP  measurements  is  smaller  than  the 
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presented  symbol  size.  The  repeatability  of  the  results  is  primarily  determined  by  the  accuracy  with  which  the 
actuator  can  be  constructed  and  adhered  to  the  surface.  We  estimate  this  as  less  than  2  mm. 

E.  Test  Conditions 

Leading  edge  separation  control  studies  have  been  performed  at  primarily  Re=Q.15  x  106  (45  m/sec)  and  1  x  106 
(62  m/sec)  with  no  flap  deflection.  Tests  are  limited  by  the  structural  fidelity  of  the  airfoil  model  which  was  not 
specifically  designed  for  high  Re ,  high  angle  of  attack  studies.  In  this  work,  the  term  “baseline”  corresponds  to  test 
cases  where  an  actuator  is  installed  on  the  surface  with  no  plasma  generation.  For  all  airfoil  measurements,  data  is 
acquired  by  establishing  a  separated  flow  baseline  condition  then  energizing  the  actuator.  For  repeated  samples  of 
different  forcing  cases  the  baseline  separated  condition  is  reestablished  between  consecutive  control  cases  to 
eliminate  hysteresis  effects. 


III.  Actuator  Characterization 


A.  Electrical  Measurements 

Figure  3a  shows  representative  discharge  pulse  voltage  and  current  traces  measured  on  a  30  cm  long  DBD 
actuator.  Figure  3b  plots  instantaneous  power  and  time-resolved  pulse  energy  coupled  to  the  load  at  these 
conditions.  During  these  measurements,  the  pulse  repetition  rate  is  fixed  at  10  Hz  to  exclude  effects  of  actuator  and 
pulse  generator  heating  on  the  coupled  pulse  energy.  The  peak  voltage  and  current  are  approximately  16  kV  and  50 
A,  respectively.  Steady  state  energy  (13  mJ/pulse)  is  realized  approximately  150  ns  after  the  beginning  of  the  pulse. 
Peak  power  of  up  to  600  kW  is  generated  due  to  the  short  pulse  width,  but  the  energy  is  less  than  1  mJ/pulse/cm  for 
all  actuator  lengths  surveyed.22  For  the  frequencies  considered  in  this  work  (<3  kHz),  this  corresponds  to  time 
average  power  of  <3  W/cm.  The  actuator  length  has  an  effect  on  these  values,  but  scaling  to  1.5  m  has  been 
accomplished  by  simply  increasing  the  DC  voltage  applied  to  the  pulse  generator.22  Note  that  the  minimum  pulse 
rise  and  fall  time  achieved  in  a  30  cm  long  actuator  is  approximately  40  nsec,  which  corresponds  to  a  wavelength  of 
X^24  m.  Therefore  at  the  present  conditions,  the  nanosecond  pulse  DBD  load  can  be  considered  as  a  lumped 
parameter  circuit  since  the  actuator  length  is  much  shorter  than  X/4.  Provided  control  authority  is  realized,  this 
slower  rise  time  is  advantageous  in  practice  because  it  permits  the  use  of  longer  actuators.  At  conditions  where  the 
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length  of  the  pulse  is  comparable  to  actuator  length,  nonuniform  potential  distribution  along  the  electrodes  becomes 
problematic.30 
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Figure  3.  Voltage  and  current  traces  (a)  and  instantaneous  power  and  energy  traces  (b)  for  a  30  cm  long  NS- 
DBD  actuator  driven  at  10  Hz. 


Figure  4  plots  the  coupled  pulse  energy  (i.e.  steady  state  energy)  against  the  pulser  operation  time  for  different 
pulse  repetition  rates  for  a  46  cm  long  actuator.  It  can  be  seen  that  the  pulse  energy  gradually  increases  with  time 
before  reaching  a  near  steady  state  value  after  approximately  10  minutes.  We  note  that  time  to  reach  the  steady  state 
is  nearly  independent  of  the  pulse  repetition  rate,  while  the  steady  state  pulse  energy  strongly  depends  on  it.  These 
results  demonstrate  that  the  energy  coupled  to  the  actuator  is  affected  significantly  by  heating.  We  believe  this  is 
primarily  due  to  heating  of  the  plasma  generator  electronics.  Details  of  this  effect  are  beyond  the  scope  of  this  paper. 

Figure  5  illustrates  the  variation  of  the  coupled  pulse  energy  at  steady  state  on  the  pulse  repetition  rate,  at  the 
conditions  of  Figure  4.  It  can  be  seen  that  increasing  the  pulse  repetition  rate  from  10  Hz  to  1  kHz  results  in  coupled 
pulse  energy  increase  by  more  than  50%,  from  8.9  mJ  to  14.8  mJ,  but  it  is  unclear  if  this  energy  increase  is  coupled 
to  the  near  surface  gas  layer.  This  is  likely  due  to  the  more  filamentary  behavior  of  the  discharge  at  high  frequency 
(see  Figure  9),  but  no  systematic  studies  have  been  undertaken  to  support  this  conjecture. 
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Figure  4.  Coupled  pulse  energy  vs.  operation  time  for  different  pulse  repetition  rates  for  a  46  cm  long 

actuator. 


Figure  5.  Coupled  pulse  energy  at  steady  state  vs.  pulse  repetition  rate  for  a  46  cm  long  actuator. 


B.  Behavior  in  Quiescent  Air 

AC-DBD  plasma  actuator’s  control  authority  is  shown  to  be  primarily  related  to  EHD  interactions  that  produce 
tangential  wall  jets  with  maximum  velocity  of  a  few  m/s.31  The  velocity  field  generated  by  the  NS-DBD  discharge  in 
still  air  has  been  examined  with  PIV.  Figure  6  presents  ensemble-averaged  velocity  profiles  for  NS-DBD  measured 
20  mm  downstream  of  the  discharge  for  repetition  frequency  of  2  kHz.  The  actuator  length  in  this  case  is  30  cm.  An 
analogous  AC-DBD  induced  velocity  profile  is  shown  for  comparison.  The  two  profiles  are  acquired  using  the  same 
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actuator  by  changing  the  input  waveform.  A  slight  wall  jet  is  generated  by  the  NS-DBD  actuator,  but  it  is  much 
weaker  than  the  AC-DBD  induced  flow.  Thus  the  EHD  effect,  while  present  for  the  NS-DBD  plasma,  is 
substantially  weaker  than  in  the  case  of  AC-DBD.  This  very  weak  EHD  effect  has  also  been  observed  using  smoke 
flow  visualization  for  a  shorter  NS-DBD  pulse  width.9  The  limitations  of  low  momentum  AC-DBD  plasma  actuators 
are  well-established.32  This  fact  with  the  results  of  Figure  6  show  that  the  EHD  effects  of  NS-DBD  are  likely  not  a 
mechanism  for  flow  control  in  the  Re  discussed  in  this  work.  Rather,  rapid  localized  heating  that  excites  natural  flow 
instabilities  is  believed  to  be  the  source  of  control  authority  similar  to  that  observed  for  LAFPAs  in  both 
experiments  and  computations.  One  of  the  characteristics  of  this  control  mechanism  is  the  generation  of 
compression  waves  due  to  rapid  local  heating  of  the  near  surface  gas  layer  by  the  plasma.  In  the  case  of  LAFPAs, 
these  waves  propagate  spherically  as  if  generated  by  point  source.15  Similar,  but  substantially  more  complex 
behavior  has  been  confirmed  for  NS-DBD  plasma  using  phase-averaged  schlieren  imaging. 


Figure  6.  Mean  velocity  profiles  for  NS  and  AC-DBD  plasma  actuators  operating  in  quiescent  air  at  2  kHz 
measured  20  mm  downstream  on  a  30  cm  DBD  load. 


Figure  7  shows  phase-averaged  schlieren  images  for  three  different  time  delays  between  the  nanosecond 
discharge  pulse  and  the  schlieren  light  source  (LED)  pulse.  Both  the  discharge  and  the  schlieren  light  source  pulse 
repetition  rate  is  1  kHz  and  the  CMOS  camera  exposure  time  156  msec  (i.e.  images  are  averaged  over  156  discharge 
pulses).  These  images  are  taken  using  an  actuator  mounted  on  a  cylinder,  looking  in  the  direction  parallel  to  the 
cylinder  axis.  The  surface  of  the  actuator  is  curved  in  the  image.  The  image  visualizes  the  density  gradient  along  the 
vertical  axis  and  the  second  derivative  of  the  density  gradient  along  the  horizontal  axis.  The  intensity  is 
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approximately  proportional  to  the  density  gradient  in  the  positive  y-direction.  Since  the  wave  generated  by  the 
discharge  pulse  is  propagating  upward,  a  lower  intensity  indicates  compression  and  a  higher  intensity  indicates 
expansion.  A  wave  system  in  Figure  7  consists  of  a  compression  wave  and  an  expansion  wave  following  it.  To 
monitor  the  wave  propagation,  a  series  of  schlieren  images  have  been  taken  for  delay  times  after  the  pulse  ranging 
from  1  psec  to  50  psec  at  1  psec  intervals.  The  compression  wave  behavior  is  relatively  complex  as  the  wavefront  is 
composed  of  both  a  quasi-planar  and  cylindrical  compression  wave,  but  appears  relatively  uniform  along  the 
viewing  axis.  This  behavior  is  qualitatively  similar  to  those  reported  by  others  in  experimental  and  numerical 
results.9,12 


Figure  7.  Phase-averaged  schlieren  images  of  a  compression  wave  generated  by  the  nanosecond  pulse 
discharge  for  several  time  delays  after  the  discharge  pulse. 

Figure  8  shows  phase-averaged  schlieren  images  taken  from  a  different  direction,  perpendicular  to  the  cylinder 
axis.  The  cylinder  is  on  the  bottom  of  the  images.  In  this  orientation  the  line-of-sight  averaged  density  gradient  is 
considerably  lower  than  in  Figure  7,  where  it  is  averaged  over  the  depth  of  field  of  approximately  10  cm.  To  combat 
this,  the  CMOS  camera  exposure  time  is  increased  and  an  additional  filter  was  used  to  remove  the  discharge 
emission.  Phase-avearged  schlieren  images  in  Figure  8  show  a  nearly  plane  wave  front  propagating  along  the 
vertical  axis,  as  well  as  multiple  circular  waves.  Similar  to  the  images  in  Figure  7,  higher  signal  intensity  indicates 
compression  and  lower  intensity  indicates  expansion.  It  can  now  be  seen  that  the  wave  front  near  the  actuator 
surface  appears  to  have  multiple  ripples  created  by  individual  circular  waves,  coalescing  into  a  superpositioned 
wave  with  an  apparently  plane  front.  It  is  likely  that  the  line-of-sight  averaging  of  these  individual  waves  creates  the 
circular  wave  front  apparent  in  schlieren  images  of  Figure  7.  Since  all  schlieren  images  shown  are  averaged  over 
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multiple  discharge  pulses,  it  appears  that  individual  waves  apparent  in  Figure  8  originate  from  relatively  stable 
filaments  which  are  repeatedly  generated  at  nearly  the  same  locations.  This  individual  circular  wave  pattern 
correlates  well  with  the  discharge  emission  intensity  pattern,  providing  additional  evidence  that  some  discharge 
filaments  remain  relatively  steady  pulse-to-pulse  and  thus  generate  repetitive  waves,  originating  at  nearly  the  same 
locations.22  This  can  also  be  visualized  using  long  exposure  photographs.  Various  pulse  repetition  rates  are  shown  in 
Figure  9.  Camera  exposure  time  is  varied  to  keep  the  number  of  discharge  pulses  during  the  exposure  approximately 
the  same.  DC  input  voltage  680  V  and  300-320  discharge  pulses  are  accumulated  in  each  image.  It  can  be  seen  that 
the  discharge  becomes  progressively  more  filamentary  as  the  pulse  repetition  rate  is  increased.  At  a  pulse  repetition 
rate  of  1 00  Hz,  the  discharge  appears  diffuse  and  nearly  uniform.  However,  as  the  repetition  rate  is  increased  above 
1  kHz,  multiple  large-scale  filaments  begin  to  form.  This  correlates  with  the  increased  energy  levels  observed  in 
Figure  4  and  5  at  higher  frequency.  The  filamentary  nature  of  NS-DBD  plasma  has  also  been  observed  in  ICCD 
images  with  nsec  exposure  time.9,22  We  postulate  that  the  spatial  location  of  the  filament  is  a  property  of  the  actuator 
construction  meaning  two  different  actuators  are  not  likely  to  develop  these  spatially  stationary  waves  at  the  same 
location.  These  strong  compression  waves  may  be  generated  if  the  rate  of  localized  heating  in  the  filaments  occurs 
on  the  time  scale  shorter  than  the  acoustic  time  scale,  Tacoustic~  d/a  -300  nsec,  where  d~ 0.1  mm  is  the  filament 
diameter  and  a- 300  m/sec  is  the  speed  of  sound.  This  rough  calculation  is  supported  by  the  employed  pulse  width 
(-100  ns). 
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Figure  8.  Phase-averaged  schlieren  images  taken  from  a  direction  perpendicular  to  the  cylinder  axis. 
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Figure  9.  Photographs  of  a  NS-DBD  at  different  pulse  repetition  rates. 


Figure  10  plots  compression  wave  front  position  above  the  cylinder  surface  against  the  time  delay  between  the 
nanosecond  discharge  pulse  and  the  schlieren  light  source  pulse.  These  data  are  obtained  from  the  schlieren  images 
of  Figure  8.  Compression  wave  front  position  is  defined  as  the  point  where  the  schlieren  signal  intensity  reaches  a 
local  maximum.  Symbol  size  in  Figure  10  indicates  the  uncertainty  in  the  wave  front  location  and  schlieren  light 
source  timing.  The  trajectory  of  the  M=1  sound  wave  at  T=285  K  (u=338  m/sec)  is  shown  as  a  straight  line  for 
reference.  One  can  see  that  a  few  mm  from  the  cylinder  surface,  the  speed  of  compression  waves  generated  by  the 
pulsed  discharge  (i.e.  the  trajectory  slope)  is  approaching  Mach  1.  We  postulate  that  the  compression  wave  strength 
near  the  discharge  may  be  greater  as  shown  by  Roupassov  et  al.9,  but  better  resolution  is  required  to  substantiate  this 
claim.  Clearly,  more  work  is  required  to  establish  a  full  characterization  for  the  amplitude  of  this  device  in  terms  of 
temperature  or  pressure  difference  across  the  wave.  In  any  case,  the  production  of  this  relatively  weak  compression 
wave  is  shown  to  have  substantial  effects  on  the  flow  field  in  Section  IV.  It  should  be  noted  that  similar  compression 
wave  behavior  does  not  exist  in  the  literature  for  AC-DBD  plasma  actuators  and  thermal  effects  from  these  appear 
insignificant.17"19 
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Figure  10.  Time-resolved  compression  wave  front  location,  inferred  from  the  phase-averaged  schlieren 
images. 


IV.  Flow  Results 


A.  Airfoil  Characterization 

The  CL>max  angle  for  the  NASA  EET  at  Re=0.15  xlO6  is  120.23  The  OSU  version  requires  a  boundary  layer  trip  to 
perform  in  a  similar  fashion.  Figure  1 1  compares  NASA  data  with  two  OSU  cases.  The  original  airfoil  with  no  trip 
installed  has  CL>max  at  10  0  with  considerably  different  dCL/da.  Application  of  a  boundary  layer  trip  to  the  leading 
edge,  in  this  case  a  passive  DBD,  corrects  the  dCL/da  discrepancy  and  produces  behavior  more  similar  to  NASA 
data.  The  stall  characteristics  are  highly  dependent  on  both  the  type  and  location  of  the  tripping  device  as  well  as  the 
wind  tunnel  flow  quality.  From  this  data,  we  conclude  that  the  separating  boundary  layer  is  nominally  laminar  (or 
possibly  transitional)  at  this  Re.  The  offset  in  CL  values  between  the  OSU  and  NASA  models  is  primarily  attributed 
to  the  lack  of  CP  resolution  in  the  OSU  case. 
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Figure  11:  Comparison  of  CL  vs.  a  curves  for  NASA23  and  OSU  EET  airfoils  for  /te=0.75  x  106. 

The  focus  of  this  work  is  on  controlling  separation  at  high  a.  This  requires  the  actuator  to  be  placed  very  near  the 
leading  edge  (x/c~ 0).  The  small  scale  model  does  not  permit  both  a  tripping  device  and  a  DBD  actuator  to  be 
installed.  The  remainder  of  this  work  employs  an  actuator  mounted  slightly  downstream  of  the  leading  edge  (6  mm 
in  arc  length)  with  plasma  formation  in  the  upstream  direction  (Figure  12).  The  motivation  for  this  is  twofold.  It 
places  the  location  of  plasma  formation  as  close  as  possible  to  the  expected  separation  line  and  eliminates  any  effect 
of  the  discontinuity  from  the  high  voltage  electrode  on  the  approaching  boundary  layer  (i.e.  the  first  surface 
discontinuity  the  flow  encounters  is  the  plasma  side  of  the  high  voltage  electrode).  Additionally,  the  dielectric  width 
is  increased  to  5.08  cm  (2  inches)  to  move  these  surface  discontinuities  farther  downstream  (Figure  12).  With  this 
arrangement,  as  the  flow  travels  from  the  stagnation  point  around  the  suction  surface  it  encounters  a  discontinuity  of 
0.09  mm  (3.5  mil)  forward  facing  step,  0.09  mm  (3.5  mil)  backward  facing  step  and  lastly  a  0.27  mm  (10.5  mil) 
backward  facing  step  (see  curved  arrow  in  Figure  12).  Actuation  (steady/unsteady  blowing  or  AC-DBD)  for  leading 
edge  separation  control  studies  is  usually  introduced  in  the  same  direction  as  the  mean  flow.  The  reverse  orientation 
is  implemented  with  the  assumption  that  the  efficacy  of  thermal  mechanisms  for  flow  control  is  not  directionally 
dependent.  It  should  be  noted  that  this  orientation  is  also  employed  by  Roupassov  et  al.9 
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Figure  12:  Photograph  of  the  NS-DBD  plasma  actuator  mounted  near  the  airfoil  leading  edge  in  reverse 
arrangement  (upstream  edge  located  6  mm  in  arc  length  from  x/c= 0). 


Any  reference  to  baseline  conditions  refers  to  the  configuration  where  an  actuator  is  mounted  to  the  leading 
edge  with  no  plasma  generation.  The  CL  characteristics  for  this  case  are  shown  in  Figure  13.  The  values  are 
calculated  by  creating  pin  holes  in  the  actuator  at  the  location  of  pressure  taps  and  integrating  CP.  This  effectively 
creates  a  short  circuit  and  prevents  plasma  generation.  The  pre-stall  performance  is  similar  to  the  trip  case  in  Figure 
11,  but  the  stall  characteristics  differ  since  this  actuator  is  not  as  effective  as  a  passive  trip.  NS-DBD  plasma 
actuation  is  employed  as  an  active  trip  at  NASA  pre-stall  a  (10-12°)  in  the  following  sections. 
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Figure  13:  Baseline  CL  vs.  a  data  for  /te=0.75  x  106. 

High  bandwidth  pressure  data  from  a  transducer  at  x/c=0.40  at  Re=0.15  x  106  are  shown  in  Figure  14.  The 
results  are  in  general  agreement  with  expectations  from  Figure  13.  PSD  magnitude  begins  to  increase  for  a=10° 
consistent  with  the  onset  of  stall.  At  higher  a,  the  PSD  magnitudes  are  further  amplified  consistent  with  full 
separation  over  the  model.  These  data  are  provided  to  show  there  are  no  coherent  dynamics  associated  with  baseline 
separation  although  there  is  broad  peak  centered  around  F*c= 2  at  a=16°. 


Figure  14:  Baseline  spectra  of  fluctuating  pressure,  cp ,  measured  at  x/c=0.40  for  Re=0J5  x  106  at  various  a. 
B.  Separation  Control 
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An  important  aspect  in  leading  edge  airfoil  separation  control  experiments  is  the  physical  mechanism  by  which 
control  is  achieved.  Some  examples  for  nominally  2D  actuation  include  accelerating  turbulent  transition  via  an 
active  trip,  exciting  spanwise  vortices  that  entrain  freestream  momentum,  and  simply  energizing  the  boundary  layer 
through  steady  blowing  allowing  it  to  negotiate  the  severe  adverse  pressure  gradient  at  the  leading  edge.  The 
actuator  characterization  data  for  NS-DBD  plasma  (Figure  6)  suggest  the  last  one  is  not  feasible  at  these  Re.  CP  and 
PIV  data  have  confirmed  the  first  two  mechanisms  can  be  accomplished  up  to  Re=  1  x  106  (62  m/s)  depending  on 
angle  of  attack  and  Re.  This  behavior  has  also  been  observed  on  a  NACA  001 5. 33 

Application  of  the  actuator  to  the  leading  edge  prevents  measurements  of  static  pressure  in  this  region  as  seen  in 
Figure  15a  for  both  a  stalled  and  attached  condition.  In  lieu  of  A CL,  the  change  in  CP  nearest  the  leading  edge  (which 
correlates  with  AQ)34  is  employed  (Figure  15a).  Frequency  sweeps  using  A CP  atx/c=0.05  as  a  metric  are  presented 
at  Re=0.15  x  106  and  1  x  106  in  Figure  15b.  Three  angles  of  attack  are  examined  (12°,  14°  and  16°)  and  each  shows  a 
different  response  to  NS-DBD  plasma  actuation.  It  should  be  noted  that  the  expected  CL,max  angle  for  this  airfoil  is 
12°  at  these  Re ,23  but  the  OSU  version  requires  a  boundary  layer  trip  to  perform  in  this  fashion  (see  Figure  11).  The 
data  corresponding  to  this  case  (a=12°)  in  Figure  15  show  no  frequency  preference  in  that  once  a  threshold  level  is 
reached,  control  authority  remains  relatively  constant.  In  the  post- stall  region  beginning  at  «=14°,  some  optimal 
frequencies  become  apparent  around  Fc+= 4.  As  a  is  further  increased,  the  preferred  frequency  increases  up  to  a 
maximum  of  around  5.6  at  a  =16°  while  the  control  authority  decreases.  The  control  authority  at  a=16°  has  been 
reduced  slightly  in  comparison  to  previous  experiments.34  This  is  expected  result  since  the  actuator  is  now  likely  aft 
of  the  optimal  location  for  control  (i.e.  separation  line)  at  a=16°.  Despite  this  reduced  control  authority  at  high  a, 
some  significant  flow  physics  are  revealed.  It  is  clear  that  the  NS-DBD  efficacy  does  not  suffer  due  to  the  increased 
Re.  In  fact,  in  the  two  lower  a  cases,  it  is  augmented.  This  is  even  more  promising  when  considering  that  dynamic 
pressure,  which  is  likely  a  more  reliable  scaling  parameter  than  velocity  or  Re  due  to  its  direct  relation  to 
momentum,  increases  in  proportion  to  Uj.  Note  also  that  Fc+  scaling  holds  quite  well  as  the  frequency  preference 
for  all  three  a  is  quite  similar.  It  is  also  clear  that  there  is  some  effect  on  the  flow  for  all  Fc+  as  the  leading  edge  CP 
value  (x/c=0.05)  has  been  improved  by  at  least  0.5  in  all  cases.  With  this  in  mind,  the  flow  physics  in  various 
frequency  regimes  is  explored  using  PIV. 
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Figure  15:  Effect  of  NS-DBD  plasma  on  suction  side  CP  at  x:/c=0.05  as  a  function  of  Fc+. 


Figure  16a  shows  the  NS-DBD  plasma  actuator  effect  on  CP  for  a=  12°.  This  coincides  with  CL>max  for  the  NASA 
airfoil  and  for  the  OSU  model  when  a  properly  designed  and  placed  trip  can  effectively  attach  the  flow.  The  forcing 
frequency  shown  corresponds  to  a  condition  where  the  threshold  frequency  has  been  reached  and  further  increases 
do  not  result  in  significant  CP  improvement.  The  baseline  flow  is  characterized  by  a  weak  adverse  pressure  gradient 
along  the  suction  surface.  This  behavior  is  apparent  in  the  ensemble-averaged  vorticity  field  in  Figure  16b  as  a 
strong  region  of  vorticity  signifies  the  boundary  between  the  freestream  and  separated  region  over  the  model. 
Control  at  Fc+= 4  fully  reattaches  the  flow  creating  a  strong  suction  peak  due  to  acceleration  at  the  leading  edge.  It  is 
also  clear  that  stagnation  point  has  moved  downstream  as  expected.  The  ensemble-averaged  vorticity  field  of  Figure 
16c  shows  stronger  levels  near  the  leading  edge  indicating  a  significant  effect  of  the  actuator  on  the  flow  in  this 
region.  Flow  separation  is  observed  near  x/c=0.8  in  both  controlled  pressure  and  vorticity  data  (Figure  16a,c). 
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Figure  16.  Time-averaged  behavior  of  C/>(a)  and  normalized  ensemble-averaged  vorticity,  Q *,  for  baseline  (b) 
and  NS-DBD  forcing  at  Fc+= 4  (c)  at  to=0.75xl06,  a= 12°. 


Normalized  swirling  strength,  Xc  *=lcic/U^  is  employed  for  vortex  identification  in  Figure  17.  Unlike  out-of¬ 
plane  vorticity  which  identifies  both  rotation  and  shear,  this  technique  ignores  regions  of  pure  shear  and  highlights 
regions  of  pure  rotation.35  The  method  is  based  on  critical  point  analysis  of  the  local  velocity  gradient  tensor  and  its 
eigenvalues.36  In  two-dimensional  form,  the  velocity  gradient  tensor  is: 


vw  = 


rdU_  dU A 
fix  dy 
dV_  dV 
v  fix  dy  j 


(1) 


where  in  this  case  the  gradient  has  been  applied  to  the  phase-averaged  velocity  field,  W  .  The  parameter  of  interest 
is  the  imaginary  component  of  the  eigenvalues  of  Equation  2,  which  is  nonzero  only  if: 
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The  imaginary  component  of  the  eigenvalues  of  Equation  2,  Xch  is  calculated  using  a  2nd  order  accurate  central 
difference  scheme  and  the  normalized  results  are  shown  in  Figure  17.  Four  phases  (AO=7i/2)  of  the  actuation  period 
are  shown.  None  of  the  images  show  any  recognizable  dynamic  content  (i.e.  coherent  structures)  although  an  effect 
of  actuation  very  near  the  FE  is  apparent.  There  is  also  weak  vortex  formation  aft  of  x/c= 0.8,  but  the  observable 
regions  of  swirl  are  mostly  stationary.  Thus,  reattachment  of  the  nominally  separated  flow  does  not  occur  through 
sustained  production  of  large  coherent  structures  in  this  case.  Instead,  the  flow  has  been  reattached  through  the  use 
of  the  actuator  as  an  active  tripping  device.  This  is  supported  by  both  NASA  and  OSU  results  (see  Figure  11). 
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Figure  17  Phase-averaged  (AO=7i/2)  normalized  swirling  strength,  for  NS-DBD  plasma  forcing  at  Fc+= 4  at 
/te=0.75xl06,  a= 12°. 
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Figure  18a  shows  the  baseline  and  controlled  CP  distributions  for  o=14°.  In  this  case,  the  baseline  is 
characterized  by  a  zero  pressure  gradient  stall  with  control  effective  in  reattaching  the  separated  flow  along  the 
suction  surface  to  approximately  70%  of  chord.  The  ensemble-averaged  normalized  vorticity  field  shows  this  effect 
as  the  region  of  strong  vorticity  bounding  the  freestream  and  separated  region  is  significantly  more  intense  than  the 
same  figure  for  the  a=12°  case  (Figure  16b).  Control  moves  the  region  of  strong  vorticity  near  the  leading  edge 
further  upstream.  This  strong  vorticity  extends  downstream  to  approximately  x/c= 0.2  after  which  a  more  diffuse 
region  is  observed  (Figure  18c).  Vorticity  for  this  controlled  case  differs  substantially  from  the  a=12°  case.  This  can 
also  be  observed  from  the  pressure  data.  The  peak  suction  value  and  stagnation  location  are  roughly  equivalent  to 
the  a=12°  case,  but  further  inspection  reveals  some  important  differences  between  the  two  pressure  distributions. 
The  value  of  CP  at  x/c=0.1  in  this  case  is  approximately  -2.6  compared  to  -2.3  in  Figure  16a.  Also,  the  average  value 
of  suction  surface  CP  fromx/c=0.7-1.0  is  -0.57  now  where  in  Figure  16a  it  is  -0.45.  Note  that  the  pressure  surface  CP 
values  are  nearly  equal  between  the  controlled  a=12°  and  14°  cases  except  at  95%  chord.  With  roughly  equivalent 
CP  at  the  leading  edge  and  along  pressure  surface,  the  decreased  pressure  along  the  suction  surface  predicts  the  CL 
value  should  be  greater  in  this  case.  However,  one  also  expects  the  drag  coefficient  to  increase  in  comparison  given 
the  suction  surface  pressure  distribution  beyond  70%  chord  which  suggests  a  thicker  downstream  wake  as  indicated 
by  the  controlled  vorticity  fields. 
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Figure  18.  Time-averaged  behavior  of  CP  (a)  and  normalized  ensemble-averaged  vorticity,  Q\  for  baseline  (b) 
and  NS-DBD  forcing  at  Fc+= 4  (c)  at  to=0.75xl06,  a= 14°. 


The  phase-averaged  data  sets  in  Figure  19  give  the  most  striking  description  of  the  controlled  flow  field. 
Coherent  structures  dominate  and  the  number  of  these  structures  is  consistent  with  forcing  at  Fc+= 4.  The 
effectiveness  of  these  vortices  for  entraining  freestream  momentum  and  reattaching  separated  flows  is  well 
established  in  the  literature.37  Four  phases  of  the  actuation  period  (AO=7i/2)  clearly  demonstrate  the  propagation  of 
these  organized  regions  along  the  airfoil  chord.  As  expected,  the  coherent  motion  begins  to  break  up  further 
downstream  (x/c>0.65).  The  dynamic  content  suggests  that,  in  contrast  to  Figure  17,  the  control  mechanism  is  not 
laminar  to  turbulent  transition,  but  rather  the  excitation  of  flow  instability  through  NS-DBD  generated  thermal 
perturbations.  It  also  explains  the  variation  in  CP  between  the  a=12°  and  14°  cases  discussed  previously  as  the  flow 
fields  are  substantially  different.  This  mechanism  is  well-established  for  separation  control  with  oscillatory 
momentum  injection,  but  to  the  authors  knowledge  this  is  the  first  evidence  of  such  behavior  using  what  we  believe 
are  primarily  thermal  perturbations.  One  also  notes  that  the  number  of  structures  visualized  over  the  chord  is 
consistent  with  forcing  at  Fc+~ 4.  This  confirms  the  belief  that  each  NS-DBD  pulse  and  its  associated  compression 
wave  generates  a  coherent  spanwise  vortex  in  these  conditions. 
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Figure  19  Phase-averaged  (AO=7i/2)  normalized  swirling  strength,  for  NS-DBD  plasma  forcing  at  Fc+= 4  at 
to=0.75xl06,  «= 14°. 
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The  frequency  preference  indicated  in  Figure  15b  for  a=16°  shows  three  relatively  distinct  forcing  regimes.  It 
should  be  noted  that  the  optimal  forcing  frequency  in  this  range  (Fc+= 5.6)  is  greater  than  those  normally  encountered 
in  such  flow  control  studies  (Fc+~1).  Some  suggest  a  more  appropriate  length  scale  may  be  the  projected  length  of 
the  airfoil  chord  (csin(a)).38  This  is  motivated  by  bluff  body  aerodynamics  where  for  example  the  cylinder  in  cross 
flow  is  scaled  by  its  diameter.  This  scaling  renders  flat  plate  airfoil  frequency  sweep  data  at  various  a  self- similar 
under  the  pretext  that  the  wake  thickness  is  approximated  by  this  simple  projection.  Indeed,  this  may  provide  scaling 
benefits  and  in  the  a=16°  case,  it  would  reduce  the  optimal  F+  to  a  value  closer  to  those  reported  in  the  literature 
(Fc+sin(a)=1.5).  However,  the  common  definition  of  F^  is  convenient  for  relating  forcing  frequency  to  the  number 
of  structures  that  exist  over  the  separated  region  thus,  we  retain  it  here.  It  has  been  suggested  that  1-4  vortices  should 
exist  over  the  controlled  region,  but  the  optimal  F^  is  dependent  on  boundary  layer  state  and  forcing  amplitude.39. 
Note  that  the  reduction  in  efficiency  often  observed  at  higher  Fc+  is  thought  to  be  related  to  the  smaller  size  of  the 
generated  vortices.  Our  data  shows  a  slightly  higher  value  (Fc+= 4-6)  than  usually  observed,  but  it  is  unlikely  that  the 
thermal  excitation  mechanism  is  the  cause  of  this  discrepancy  since  the  flow  physics  (i.e.  production  of  spanwise 
vortices)  are  consistent  with  literature.  This  is  supported  by  recent  experiments  on  a  NACA  0015  with  NS-DBD  that 
indicate  optimal  F^  near  2. 33  Certainly,  F^-l  provides  guidance  for  selecting  forcing  frequencies,  but  is  by  no  means 
a  universal  optimum. 

Returning  to  the  CP  and  PIV  data,  low  frequency  forcing  (Fc+=0.6)  has  an  effect  on  the  naturally  stalled  flow  as 
seen  in  the  CP  curves,  but  in  this  case  the  minimum  value  is  not  observed  nearest  the  leading  edge  (Figure  20a). 
Rather  it  appears  further  downstream  near  x/c=0.15.  The  ensemble-averaged  vorticity  certainly  shows  this  change  in 
behavior  as  the  previously  apparent  shear  layer  is  diffused  considerably  (Figure  20b, c).  As  before,  the  dynamic 
content  of  the  flow  is  most  compelling  (Figure  21).  In  this  low  frequency  case,  regions  of  rotation  are  massively 
amplified  in  scale  as  they  move  over  the  model  surface.  Note  that  these  structures  do  not  follow  the  airfoil  surface  as 
closely  as  those  shown  in  Figure  19.  Instead,  they  are  swept  downstream  by  freestream  momentum  rather  quickly 
once  they  move  away  from  the  leading  edge.  The  magnitude  of  velocity  fluctuations  associated  with  this  motion  is 
approximately  twice  as  large  as  similar  data  associated  with  Figure  19.34  These  massive  coherent  structures  have 
been  observed  in  all  the  low  frequency  forcing  cases  in  Figure  15b.  The  vortex  shedding  associated  with  this  low 
frequency  forcing  has  also  has  strong  audible  signature  which  is  associated  with  the  vortex  and  not  the  plasma  itself. 
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Figure  20.  Time-averaged  behavior  of  CP  (a)  and  normalized  ensemble-averaged  vorticity,  Q\  for  baseline  (b) 


and  NS-DBD  forcing  at  Fc+= 0.6  (c)  at  /?e=0.75xl0(',  a= 16° 
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Figure  21  Phase-averaged  (AO=7i/2)  normalized  swirling  strength,  for  NS-DBD  plasma  forcing  at  Fc+=0.6 
at  /te=0.75x!06,  a= 16°. 
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Optimal  forcing  (Fc+= 5.6)  at  a=  16°  produces  reattachment  over  approximately  half  the  model  chord  (Figure 
22a).  This  can  be  identified  by  an  adverse  pressure  gradient  that  exists  up  to  x/c~ 0.6.  The  control  authority  in  this 
case  is  diminished  compared  to  optimal  forcing  at  a=14°  (Figure  18a)  as  the  suction  peak  only  reaches  CP=- 3.  The 
pressure  side  CP  is  now  slightly  greater  than  observed  in  Figure  18a,  but  the  value  on  the  suction  surface  is 
considerably  less  with  an  average  of  -0.8  after  x/c=0A.  This  suggests  that  the  CL  is  still  quite  high,  but  CD  is 
becoming  problematic  due  to  a  thicker  wake.  The  ensemble-averaged  vorticity  further  elucidates  the  effects  of 
control.  Forcing  substantially  changes  the  flow  field  as  the  intense  shear  layer  over  the  entire  model  is  eliminated 
(Figure  22b, c).  However,  evidence  of  thicker  wake  in  comparison  to  Figure  18c  is  apparent  as  indicated  from 
pressure  data.  The  reduced  control  authority  in  this  case  is  likely  due  to  the  movement  of  the  separation  line 
upstream  resulting  in  a  less  than  optimized  position  for  the  plasma  actuator. 


Figure  22.  Time-averaged  behavior  of  CP  (a)  and  normalized  ensemble-averaged  vorticity,  Q\  for  baseline  (b) 
and  NS-DBD  forcing  at  Fc+= 5.6  (c)  at  /?e=0.75xl0(',  a= 16°. 

The  phase-averaged  velocity  fluctuations  again  reveal  clear  dynamic  content  with  scales  indicative  of 
forcing  at  Fc+~ 5.6.  The  coherent  structures  do  not  propagate  as  far  downstream  as  those  at  a=14°  before 
disintegrating.  However,  their  organization  near  the  leading  edge  is  apparent  in  the  phase-averaged  swirling  strength 
which  shows  clear  vortex  formation  and  convection  near  the  leading  edge  (x/c<0.2)  (Figure  23). 
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Figure  23  Phase-averaged  (AO=7i/2)  normalized  swirling  strength,  for  NS-DBD  plasma  forcing  at  Fc+= 5.6 
at/te=0.75x!06,  a= 16°. 
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Finally  the  effects  of  high  frequency  forcing  (Fc+=  11.3)  are  investigated.  Although  the  leading  edge  CP  value  is 
similar  to  the  low  frequency  forcing  case  (Figure  20a),  the  form  of  the  CP  curve  just  downstream  of  the  leading  edge 
suggests  a  very  different  flow  field  (Figure  24a).  As  in  all  cases,  the  ensemble-averaged  vorticity  is  significantly 
altered  by  control  (Figure  24b, c),  but  merely  breaking  up  the  intense  shear  layer  over  the  airfoil  is  not  sufficient  to 
substantially  reduce  separation.  The  dynamic  features  are  now  much  less  apparent,  but  still  exist  near  the  leading 
edge  (Figure  25)  highlighting  the  robustness  of  the  swirling  strength  metric.  Despite  the  weak  oscillations,  distinct 
vortex  formation  and  convection  is  identified  at  x/c<0.1  at  a  wavelength  approximately  half  of  that  observed  in 
Figure  23. 


Figure  24.  Time-averaged  behavior  of  CP  (a)  and  normalized  ensemble-averaged  vorticity,  £2\  for  baseline  (b) 
and  NS-DBD  forcing  at  Fc+=  11.3  (c)  at  to=0.75xl06,  a= 16°. 
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Figure  25  Phase-averaged  (AO=7i/2)  normalized  swirling  strength,  2a*,  for  NS-DBD  plasma  forcing  at 
Fc+=11.3  at  Fe=0.75xl06,  «= 16°. 
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The  preceding  results  for  7te=0.75xl06  have  also  been  acquired  at  ite=lxl06.  As  the  results  in  Figure  15b 
indicate,  the  flow  fields  are  quite  similar  in  terms  of  control  authority,  ensemble-averaged  and  phase-averaged 
contents.  The  efficacy  at  which  the  NS-DBD  plasma  actuator  excites  coherent  structures  at  the  airfoil  leading  edge  is 
particularly  appealing.  Although  optimal  frequencies  exist  for  control,  it  seems  that  the  NS-DBD  possesses 
sufficient  authority  to  manipulate  flow  instabilities  and  create  structures  at  any  frequency  in  the  range  of  flow 
conditions  surveyed.  This  is  consistent  with  the  high  amplitude  characteristics  of  LAFPAs. 13,15 

The  CL  curve  for  the  controlled  cases  is  shown  in  Figure  26  as  a  final  metric  for  demonstrating  NS-DBD  plasma 
actuator’s  capability.  The  baseline  begins  to  stall  around  a=10°  at  Re=0.75xl06  (45  m/s)  and  the  NS-DBD  plasma 
extends  the  stall  angle  up  to  16°.  The  dimensionless  excitation  frequency  at  a=10°  is  4.  CL  values  for  a=12°,  14°  and 
16°  have  been  calculated  from  the  data  in  Figures  16,  18  and  22,  respectively.  Recall  that  CP  values  near  the  leading 
edge  are  not  measurable  due  to  covering  of  the  pressure  taps  by  the  actuator.  Consequently,  the  CL  values  have  been 
estimated  using  the  measured  CP  values  and  linearly  extrapolating  the  NS-DBD  controlled  data  at  the  leading  edge 
based  on  pre-stall  values  for  a  passive  (no  plasma)  actuator  modified  to  allow  CP  measurements.  This  methodology 
is  similar  to  other  simplified  CL  estimates,40  but  also  employs  the  measured  data  away  from  the  airfoil  leading  edge. 
The  error  bars  on  the  figure  are  intentionally  large  due  to  the  uncertainty  introduced  by  this  extrapolation.  As 
suggested  by  previous  data,  the  NS-DBD  is  capable  of  extending  the  stall  angle  to  16°,  but  we  make  no  claims  on 
the  accuracy  of  the  CL  values.  Rather,  the  key  result  here  is  that  in  terms  of  stall  angle  only,  the  NS-DBD  rivals  the 
performance  of  a  traditional  passive  control  device  (i.e.  leading  edge  droop).23 


Figure  26:  Effect  of  active  control  with  NS-DBD  plasma  actuation  on  CL  for  the  OSU  airfoil. 
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V.  Discussion  of  Results 


With  details  of  the  flow  physics  now  presented,  it  is  appropriate  to  return  to  the  actuator  characterization  results 
for  discussion  of  the  connection  between  them.  In  Figure  4,  a  transient  effect  believed  to  be  associated  with  heating 
of  the  pulser  electronics  is  indicated.  This  transient  phenomenon  persists  for  approximately  10  minutes  before 
reaching  steady  state.  During  this  time,  the  coupled  pulse  energy  is  increasing.  The  data  presented  in  Section  IV 
have  been  acquired  during  the  transient  portion  of  the  curve.  For  example,  a  set  of  PIV  data  requires  -75  seconds  for 
acquisition  (750  samples  at  10  Hz).  No  obvious  change  in  control  authority  is  observed  during  a  given  experiment 
while  monitoring  CP  data  at  1  Hz  suggesting  this  transient  has  little  effect  on  the  flow  control  results,  but  these  have 
not  yet  been  extended  past  the  1 0  min  interval.  The  degradation  of  an  actuator  will  eventually  play  a  role  as  run  time 
and  run  frequency  increase.  A  similar  situation  is  found  in  Figure  5  which  shows  that  the  coupled  energy  increases 
with  pulse  repetition  frequency.  We  postulate  that  this  is  due  to  the  more  filamentary  behavior  shown  in  Figure  9.  At 
this  time  it  is  unclear  if  the  increase  in  coupled  energy  results  in  increased  heating  of  the  near  surface  gas  layer  (and 
ultimately  the  strength  of  the  compression  wave)  or  merely  heating  of  the  actuator  material.  The  data  in  Figure  15b 
do  not  show  a  similar  correlation  as  that  of  Figure  9  as  higher  frequency  does  not  result  in  more  control  authority 
and  Fc+  scaling  is  dominant.  This  is  most  strikingly  depicted  in  Figure  21  where  a  relatively  low  frequency  (-100 
Hz)  (i.e.  low  pulse  energy)  has  a  dramatic  effect  on  the  flow  field.  A  calculation  of  the  efficiency  and  amplitude 
characteristics  of  the  device  requires  detailed  near  surface  measurement  of  temperature  rise  and/or  compression 
wave  strength  with  and  without  the  effects  of  freestream  and  efforts  to  this  aim  are  ongoing.22  In  any  case,  it  is  clear 
that  the  velocity  of  the  generated  compression  wave  becomes  approximately  sonic  a  few  mm  from  the  discharge 
(Figure  10).  It  is  perhaps  surprising  that  these  weak  waves  generate  the  high  amplitude  characteristics  observed  in 
Figure  20-25,  but  we  note  that  similar  wave  speeds  have  been  found  in  the  case  of  LAFPAs15  and  higher  near 
surface  wave  velocities  for  NS-DBD  plasma  have  been  suggested.9,12 

The  schlieren  images  of  the  Figure  8  show  fairly  complex  3D  compression  wave  behavior,  yet  the  flow  physics 

reveal  coherent  spanwise  vortices.  These  vortices  are  consistently  generated  for  physically  different  actuators 

applied  to  the  airfoil  leading  edge.  Thus,  the  flat  wave  front  found  in  the  schlieren  images  in  Figure  8  appears 

sufficient  for  producing  spanwise  vortices  here.  The  3D  structure  of  the  discharge  and  its  effect  on  the  flow  are 

intended  for  future  study.  These  results  show  compelling  evidence  for  the  use  of  repetitive  nanosecond  pulses 

applied  to  a  DBD  load  for  separation  control.  This  generally  supports  the  findings  of  Roupassov  et  al.9  at  takeoff  and 
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landing  velocities  associated  with  typical  transport  aircraft  although  some  distinctions  can  be  made.  The  cited  work 
employed  shorter  pulse  width  (-5-70  ns)  and  higher  pulse  voltage  (25  kV)  to  guarantee  the  necessary  heating  in  the 
near  surface  gas  layer.  The  pulse  width  (-100  ns)  and  pulse  voltage  (-15  kV)  employed  in  our  work  is  sufficient  to 
generate  substantial  control  authority  at  the  tested  freestream  conditions.  It  is  probably  no  coincidence  that  the  pulse 
energy  between  these  two  studies  is  of  the  same  order  of  magnitude  (-0.5  mJ/cm). 


VI.  Conclusions 

The  efficacy  of  repetitive  nanosecond  pulse  driven  dielectric  barrier  discharge  (NS-DBD)  plasma  actuators  for 
controlling  flow  separation  from  an  airfoil  leading  edge  has  been  demonstrated  up  to  Re=  lxlO6  (62  m/s).  Results 
generally  support  the  findings  of  Roupassov  et  al.9  but  the  employed  pulse  parameters  such  as  pulse  width  and  peak 
voltage  are  less  restrictive.  The  nanosecond  pulse  plasma  generator  used  in  this  work  has  been  custom  designed  and 
built  in-house.  The  device  has  been  characterized  in  quiescent  air  at  atmospheric  pressure.  Current,  voltage, 
instantaneous  power  and  coupled  pulse  energy  have  been  measured  for  different  pulse  peak  voltages,  pulse 
repetition  rates,  and  actuator  lengths.  The  results  show  that  the  actuator  can  be  scaled  to  length  of  at  least  1.5  m  by 
simply  increasing  the  DC  supply  voltage.22  The  coupled  energy  of  the  plasma  is  less  than  1  mJ/pulse/cm  in  all  cases 
(power  <3  W/cm)  which  is  sufficient  for  producing  control  authority  in  this  work.  This  energy  is  dependent  on  both 
the  run  time  and  frequency  of  the  discharge,  but  control  authority  does  not  show  obvious  dependence  on  coupled 
energy.  Rather  it  obeys  T7'  scaling,  presumably  because  sufficient  amplitude  is  produced  even  for  low  pulse  energy. 

The  NS-DBD  transfers  substantially  less  momentum  to  the  neutral  air  molecules  than  the  more  widely  studied 
AC-DBD  plasma  actuator.  The  momentum  transfer  is  quite  weak  and  is  not  expected  to  have  a  significant  effect  in 
the  flow  in  the  conditions  tested.  Instead,  the  control  authority  generated  by  NS-DBD  plasma  is  believed  to  stem 
from  rapid  localized  heating  of  the  near  surface  gas  layer  by  the  plasma.  This  gives  rise  to  a  complex  pattern  of 
quasi-planar  and  spherical  compression  waves  in  still  air.  Phase-averaged  schlieren  images  suggest  that  some  of 
these  compression  waves  are  generated  by  discharge  filaments  that  remain  fairly  reproducible  pulse-to-pulse. 
Despite  the  3D  behavior  observed  in  quiescent  air,  a  nominally  2D  effect  is  observed  on  the  airfoil  leading  edge  as 
coherent  spanwise  vortices  are  developed.  At  nominally  pre-stall  angles  of  attack  for  high  Re  (a=  12°),  the  NS-DBD 
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functions  as  an  active  trip  while  at  post-stall  it  is  found  to  manipulate  flow  instabilities  to  generate  coherent 
spanwise  vortices  across  the  frequency  range  surveyed.  The  latter  is  widely  recognized  for  controlling  separation  in 
a  variety  of  flow  systems.  The  device  appears  to  possess  high  amplitude  characteristics  as  structures  are  generated  at 
all  surveyed  frequencies,  but  Fc+  values  in  the  range  4-6  are  most  effective  for  controlling  flow  separation.  The 
control  authority  ofNS-DBD  does  not  suffer  from  Re  effects  between  7te=0.75xl06  (45  m/s)  and  Re=  lxlO6  (62  m/s) 
and  follows  Fc+  scaling  as  expected.  At  this  time,  the  bandwidth  of  the  device  is  limited  to  <3  kHz  due  to  heating  of 
the  power  supply,  but  there  is  no  lower  limit  on  pulse  repetition  frequency.  Efforts  to  extend  the  operating  frequency 
are  ongoing.  It  should  also  be  emphasized  that  burst  mode  of  operation  has  not  been  employed  and  that  each  NS- 
DBD  pulse  is  found  to  generate  these  large  scale  structures. 

This  work  presents  further  evidence  that  the  thermal  mechanism  observed  in  LAFPAs13,15  and  also  presented  by 
Roupassov  et  al.9  can  be  employed  for  separation  control  with  a  DBD  load.  This  may  permit  the  use  of  DBD  plasma 
actuators  for  high  speed  flow  control  applications  (Re>lxl06,  f/^100  m/s)  although  the  freestream  conditions 
presented  here  must  be  extended  to  verify  this  hypothesis.  Substantial  work  is  also  required  to  examine  the  effect  of 
boundary  layer  state  and  thickness  on  NS-DBD  control  authority. 
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